Dissection of Cellular Communication between Human Primary Osteoblasts and Bone Marrow Mesenchymal Stem Cells in Osteoarthritis at Single-Cell Resolution

Background and Objectives@#Osteoblasts are derived from bone marrow mesenchymal stem cells (BMMSCs) and playimportant role in bone remodeling. While our previous studies have investigated the cell subtypes and heterogeneity in osteoblasts and BMMSCs separately, cell-to-cell communications between osteoblasts and BMMSCs in vivo in humans have not been characterized. The aim of this study was to investigate the cellular communication between human primary osteoblasts and bone marrow mesenchymal stem cells. @*Methods@#and Results: To investigate the cell-to-cell communications between osteoblasts and BMMSCs and identifynew cell subtypes, we performed a systematic integration analysis with our single-cell RNA sequencing (scRNA-seq) transcriptomes data from BMMSCs and osteoblasts. We successfully identified a novel preosteoblasts subtype which highly expressed ATF3, CCL2, CXCL2 and IRF1. Biological functional annotations of the transcriptomes suggested that the novel preosteoblasts subtype may inhibit osteoblasts differentiation, maintain cells to a less differentiated status and recruit osteoclasts. Ligand-receptor interaction analysis showed strong interaction between mature osteoblasts and BMMSCs. Meanwhile, we found FZD1 was highly expressed in BMMSCs of osteogenic differentiation direction. WIF1 and SFRP4, which were highly expressed in mature osteoblasts were reported to inhibit osteogenic differentiation. We speculated that WIF1 and sFRP4 expressed in mature osteoblasts inhibited the binding of FZD1 to Wnt ligand in BMMSCs, thereby further inhibiting osteogenic differentiation of BMMSCs. @*Conclusions@#Our study provided a more systematic and comprehensive understanding of the heterogeneity of osteogenic cells. At the single cell level, this study provided insights into the cell-to-cell communications between BMMSCs and osteoblasts and mature osteoblasts may mediate negative feedback regulation of osteogenesis process.

Mass spectrometry based proteomics profiling of human monocytes

Human monocyte is an important cell type which is involved in various complex human diseases. To better understand the biology of human monocytes and facilitate further studies, we developed the first comprehensive proteome knowledge base specifically for human monocytes by integrating both in vivo and in vitro datasets. The top 2000 expressed genes from in vitro datasets and 779 genes from in vivo experiments were integrated into this study. Altogether, a total of 2237 unique monocyte-expressed genes were cataloged. Biological functions of these monocyte-expressed genes were annotated and classified via Gene Ontology (GO) analysis. Furthermore, by extracting the overlapped genes from in vivo and in vitro datasets, a core gene list including 541 unique genes was generated. Based on the core gene list, further gene-disease associations, pathway and network analyses were performed. Data analyses based on multiple bioinformatics tools produced a large body of biologically meaningful information, and revealed a number of genes such as SAMHD1, G6PD, GPD2 and ENO1, which have been reported to be related to immune response, blood biology, bone remodeling, and cancer respectively. As a unique resource, this study can serve as a reference map for future in-depth research on monocytes biology and monocyte-involved human diseases.

Replication of Caucasian Loci Associated with Osteoporosis-related Traits in East Asians

BACKGROUND: Most reported genome-wide association studies (GWAS) seeking to identify the loci of osteoporosis-related traits have involved Caucasian populations. We aimed to identify the single nucleotide polymorphisms (SNPs) of osteoporosis-related traits among East Asian populations from the bone mineral density (BMD)-related loci of an earlier GWAS meta-analysis. METHODS: A total of 95 SNPs, identified at the discovery stage of the largest GWAS meta-analysis of BMD, were tested to determine associations with osteoporosis-related traits (BMD, osteoporosis, or fracture) in Korean subjects (n=1,269). The identified SNPs of osteoporosis-related traits in Korean subjects were included in the replication analysis using Chinese (n=2,327) and Japanese (n=768) cohorts. RESULTS: A total of 17 SNPs were associated with low BMD in Korean subjects. Specifically, 9, 6, 9, and 5 SNPs were associated with the presence of osteoporosis, non-vertebral fractures, vertebral fractures, and any fracture, respectively. Collectively, 35 of the 95 SNPs (36.8%) were associated with one or more osteoporosis-related trait in Korean subjects. Of the 35 SNPs, 19 SNPs (54.3%) were also associated with one or more osteoporosis-related traits in East Asian populations. Twelve SNPs were associated with low BMD in the Chinese and Japanese cohorts. Specifically, 3, 4, and 2 SNPs were associated with the presence of hip fractures, vertebral fractures, and any fracture, respectively. CONCLUSIONS: Our results identified the common SNPs of osteoporosis-related traits in both Caucasian and East Asian populations. These SNPs should be further investigated to assess whether they are true genetic markers of osteoporosis.

Genome wide association study: searching for genes underlying body mass index in the Chinese / 生物医学与环境科学（英文）

<p><b>OBJECTIVE</b>Obesity is becoming a worldwide health problem. The genome wide association (GWA) study particularly for body mass index (BMI) has not been successfully conducted in the Chinese. In order to identify novel genes for BMI variation in the Chinese, an initial GWA study and a follow up replication study were performed.</p><p><b>METHODS</b>Affymetrix 500K SNPs were genotyped for initial GWA of 597 Northern Chinese. After quality control, 281,533 SNPs were included in the association analysis. Three SNPs were genotyped in a Southern Chinese replication sample containing 2 955 Chinese Han subjects. Association analyses were performed by Plink software.</p><p><b>RESULTS</b>Eight SNPs were significantly associated with BMI variation after false discovery rate (FDR) correction (P=5.45×10⁻⁷-7.26×10⁻⁶, FDR q=0.033-0.048). Two adjacent SNPs (rs4432245 & rs711906) in the eukaryotic translation initiation factor 2 alpha kinase 4 (EIF2AK4) gene were significantly associated with BMI (P=6.38×10⁻⁶ & 4.39×10⁻⁶, FDR q=0.048). In the follow-up replication study, we confirmed the associations between BMI and rs4432245, rs711906 in the EIF2AKE gene (P=0.03 & 0.01, respectively).</p><p><b>CONCLUSION</b>Our study suggests novel mechanisms for BMI, where EIF2AK4 has exerted a profound effect on the synthesis and storage of triglycerides and may impact on overall energy homeostasis associated with obesity. The minor allele frequencies for the two SNPs in the EIF2AK4 gene have marked ethnic differences between Caucasians and the Chinese. The association of the EIF2AK4 gene with BMI is suggested to be 'ethnic specific' in the Chinese.</p>

Genome-wide Association Studies for Osteoporosis: A 2013 Update

In the past few years, the bone field has witnessed great advances in genome-wide association studies (GWASs) of osteoporosis, with a number of promising genes identified. In particular, meta-analysis of GWASs, aimed at increasing the power of studies by combining the results from different study populations, have led to the identification of novel associations that would not otherwise have been identified in individual GWASs. Recently, the first whole genome sequencing study for osteoporosis and fractures was published, reporting a novel rare nonsense mutation. This review summarizes the important and representative findings published by December 2013. Comments are made on the notable findings and representative studies for their potential influence and implications on our present understanding of the genetics of osteoporosis. Potential limitations of GWASs and their meta-analyses are evaluated, with an emphasis on understanding the reasons for inconsistent results between different studies and clarification of misinterpretation of GWAS meta-analysis results. Implications and challenges of GWAS are also discussed, including the need for multi- and inter-disciplinary studies.

ANKRD7 and CYTL1 are novel risk genes for alcohol drinking behavior / 中华医学杂志(英文版)

<p><b>BACKGROUND</b>Alcohol dependence (AD) is a complex disorder characterized by impaired control over drinking. It is determined by both genetic and environmental factors. The recent approach of genome-wide association study (GWAS) is a powerful tool for identifying complex disease-associated susceptibility alleles, however, a few GWASs have been conducted for AD, and their results are largely inconsistent. The present study aimed to screen the loci associated with alcohol-related phenotypes using GWAS technology.</p><p><b>METHODS</b>A genome-wide association study with the behavior of regular alcohol drinking and alcohol consumption was performed to identify susceptibility genes associated with AD, using the Affymetrix 500K SNP array in an initial sample consisting of 904 unrelated Caucasian subjects. Then, the initial results in GWAS were replicated in three independent samples: 1972 Caucasians in 593 nuclear families, 761 unrelated Caucasian subjects, and 2955 unrelated Chinese Hans.</p><p><b>RESULTS</b>Several genes were associated with the alcohol-related phenotypes at the genome-wide significance level, with the ankyrin repeat domain 7 gene (ANKRD7) showing the strongest statistical evidence for regular alcohol drinking and suggestive statistical evidence for alcohol consumption. In addition, certain haplotypes within the ANKRD7 and cytokine-like1 (CYTL1) genes were significantly associated with regular drinking behavior, such as one ANKRD7 block composed of the SNPs rs6466686-rs4295599-rs12531086 (P = 6.51 × 10(-8)). The association of alcohol consumption was successfully replicated with rs4295599 in ANKRD7 gene in independent Caucasian nuclear families and independent unrelated Chinese Hans, and with rs16836497 in CYTL1 gene in independent unrelated Caucasians. Meta-analyses based on both the GWAS and replication samples further supported the observed significant associations between the ANKRD7 or CYTL1 gene and alcohol consumption.</p><p><b>CONCLUSION</b>The evidence suggests that ANKRD7 and CYTL1 genes may play an important role in the variance in AD risk.</p>

Contribution of genotype and ethnicity to bone mineral density variation in Caucasians and Chinese: a test for five candidate genes for bone mass / 中华医学杂志(英文版)

<p><b>BACKGROUND</b>Ethnicity is shown to be one of important factors affecting bone mineral density (BMD). The present study was performed to compare the association of six markers for five candidate genes with BMD variation in two populations of different ethnicity, Caucasian and Chinese, and the contribution of genotype and ethnicity to this variation in the populations.</p><p><b>METHODS</b>The studied restriction fragment length polymorphisms were BsaH I of the calcium-sensing receptor gene, SacI of the alpha2HS-glycoprotein (AHSG) gene, PvuII and XbaI of the oestrogen receptor alpha gene, ApaI of the vitamin D receptor (VDR) gene and BstBI of the parathyroid hormone gene. The association of these markers with BMD was analysed by one-way and two-way ANOVA with adjustment for covariates.</p><p><b>RESULTS</b>Two polymorphisms, AHSG-SacI and VDR-ApaI, showed no association with BMD, while the others were associated with BMD variation at some skeletal sites in either males or females. The polymorphisms indicated clear distinctions between the associations depending on ethnicity, gender and skeletal site. Similar patterns were observed in their contribution to the total population BMD variation. Ethnicity appears to have a larger effect on the total population BMD variation in females than in males. It may account, on the average, for about 2% total population BMD variation at the spine of females and about 1% at the hip of males and females.</p><p><b>CONCLUSION</b>The results of the present study suggest that significant interethnic differentiation at some loci may contribute to the significant interethnic difference in BMD. However, this contribution apparently is not large.</p>

Preliminary linkage analysis of a Chinese family with benign familial infantile convulsion / 中华儿科杂志

<p><b>OBJECTIVE</b>Benign familial infantile convulsions (BFIC) is a recently recognized autosomal dominant inherited disorder. This epileptic syndrome typically begins between 3 and 12 months of age with clusters of partial seizures in most cases and carries a good prognosis. So far, three loci have been linked to chromosome 19q12.1-13.1, chromosome 2q24 and chromosome 16p12-q12. The authors performed linkage analysis on this pedigree.</p><p><b>METHODS</b>A four-generation Chinese family was investigated. The total number of members was 32 in this family and two neurologists in Xiangya Hospital gave systemic physical examinations and interictal neurological examinations to nineteen members of this family. Venous blood samples were taken for genetic analysis. DNA was extracted from peripheral blood leukocytes using phenol-chloroform method. Seventeen microsatellite markers spanning the critical regions on chromosomes 19q12-13.1, 2q24, and 16p12-q12 were genotyped. These markers included D19S49, D19S250, D19S414, D19S416 and D19S245 for the 19q region, D2S2380, D2S399, D2S111, D2S2195, D2S2330 and D2S2345 for the 2q region, D16S401, D16S3131, D16S3093, D16S517, D16S3120 and D16S415 for the 16p-q region. The DNA from each sample was amplified for the 17 markers. After polymerase chain reactions (PCR), PCR products of chromosome 19 with markers D19S49, D19S250, D19S414, D19S416 and D19S245 were subjected to electrophoresis on 8% denatured polyacrylamide gel for at least 2 hours and 20 minutes. Then the length of the PCR products was judged in the Strategene Eagle Eye II automated gel image analyzer. For the markers from chromosome 2 and 16, PCR products were scanned at ABI 377 autosequencer. The data of PCR products were analyzed using the software Genescan v3.1, Genetyper v2.1 (Applied Biosystem, CA. USA) and GenoDB v1.0. After Mendelian checking, the eligible genotyping data were used for linkage analysis. LOD scores were calculated by using MLINK program of LINKAGE v5.1, under an assumption of autosomal dominant inheritance and the estimated penetrance was 0.9. The allele frequencies of each marker were assumed to be equal and the disease-allele frequencies were designated to be 1/10,000. The LOD scores were calculated at combination rate (theta) 0.0, 0.1, 0.2, 0.3, and 0.4.</p><p><b>RESULTS</b>Among the 17 selected microsatellite markers, which cover the previously reported regions, seven markers' data (D16S3131, D16S517, D16S3120, D16S3093, D2S2380, D19S250 and D19S414) were omitted due to failed genotyping, low genetic heterogeneity, or failure to pass Mendelian checking. Omission of these markers was to ensure the reliability of our raw data. The two-point LOD scores were below zero for all the markers and the maximum LOD scores at theta = 0.0 were less than -2 for markers D19S49, D19S416, D19S245, D16S401, D16S415, D2S399, D2S111, D2S2195, D2S2330 and D2S2345. Thus, the linkage result showed no evidence that the disease locus is linked to any of these selected markers, which excludes the previously reported candidate regions found in other ethnic families.</p><p><b>CONCLUSION</b>There is no evidence that this Chinese family was linked to one of the following loci: 19q12.1-13.1, 16p12-q12 and 2q24. The results indicated that BFIC showed genetic heterogeneity and the Chinese BFIC families might be mapped on another new locus.</p>

Fine mapping of susceptibility genes by Lewontin's linkage disequilibrium measure with application to Alzheimer's disease / 中华医学杂志(英文版)

<p><b>OBJECTIVES</b>To formulate an equation for fine mapping of disease loci under complex conditions and determine the marker-disease distance in a specific case using this equation.</p><p><b>METHODS</b>Lewontin's linkage disequilibrium (LD) measure D' was used to formulate an equation for mapping disease genes in the presence of phenocopies, locus heterogeneity, gene-gene and gene-environment interactions, incomplete penetrance, uncertain liability and threshold, incomplete initial LD, natural selection, recurrent mutation, high disease allele frequency and unknown mode of inheritance. This equation was then used to determine the distance between a marker ( epsilon 4 within the apolipoprotein E gene, APOE) and Alzheimer's disease (AD) loci using published data.</p><p><b>RESULTS</b>An equation was formulated for mapping disease genes under the above conditions.If these conditions are present but ignored, then recombination fraction theta between marker and disease loci will be either overestimated or estimated with little bias. Therefore, an upper limit of theta can be obtained. AD has been found to be associated with the marker allele epsilon 4 in Africans, Asians, and Caucasians. This suggests that the AD- epsilon 4 allelic LD predates the divergence of peoples occurring 100 000 years ago. With the age of AD- epsilon 4 allelic LD so estimated, the maximal distance was calculated to be 23.2 kb (mean 5.8 kb).</p><p><b>CONCLUSIONS</b>(1) A method is developed for LD mapping of susceptibility genes. (2) A mutation within the APOE gene itself, among others, is responsible for the susceptibility to AD, which is supported by recent evidence from studies using transgenic mice.</p>